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^TT* 1 ■ Abstract 

O 

We apply Christodoulou's framework, developed to study the Einstein-scalar field equations in 

spherical symmetry, to the linear wave equation in de Sitter spacetime, as a first step towards the 
Einstein-scalar field equations with positive cosmological constant. We obtain an integro-differential 
evolution equation which we solve by taking initial data on a null cone. As a corollary we obtain 
elementary derivations of expected properties of linear waves in de Sitter spacetime: boundedness 
w~, . in terms of (characteristic) initial data, and a Price law establishing uniform exponential decay, in 
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Bondi time, to a constant. 



^ 1 Introduction 

The study of the linear wave equation 



□ 9 = O (1) 



p^j . on fixed backgrounds (M, g) is a stepping stone to the analysis of the nonlinearities of gravitation. In 

J> " this paper we apply Christodoulou's framework, developed in 3 , to spherically symmetric solutions 

of |T|) on a de Sitter background, as a prerequisite to the study of the coupled Einstein-scalar field 
equations with positive cosmological constant in spherical symmetry (which will be pursued elsewhere). 
If the cosmological constant vanishes then the uncoupled problem ([TJ is trivial]], a fact that Christodoulou 
explored in [3] to solve the coupled case for suitably small initial data. For positive cosmological constant, 
however, the uncoupled case is more complicated^, and it is essential to understand it thoroughly in order 
to ascertain how much freedom is there when perturbing it to the nonlinear case, as well as to determine 
which decays to expect and which function spaces to use. 

Following Christodoulou, we turn ((TJ into an integro-differential evolution equation, which we solve 
by taking initial data on a null cone. This step, which is trivial in the case of vanishing cosmological 
constant, turns out to be quite subtle for positive cosmological constant. As a corollary we obtain 
elementary derivations of expected properties of linear waves in de Sitter spacetime: boundedness in 
terms of (characteristic) initial data, and uniform exponential decay, in Bondi time, to a constant (from 
which exponential decay to a constant in the usual static time coordinate easily follows; such boundedness 
and decay results may be seen, respectively, as analogues of the Kay-Wald theorem [5] and of a Price 
law jTUJ |3] , both originally formulated for Cauchy data in a Schwarzschild background) . Although widely 
expected from the stability of de Sitter spacetime [7J, we are unaware of a written proof of uniform 
exponential decay. Similarly, the bound that we obtain for the solution in terms of the C° norm of the 
(characteristic) initial data is, to the best of our knowledge, original (notice that in particular the bound 
([7]) involves no "loss of derivatives"). Another novel aspect of our work is the fact that our results apply 
to a domain containing the cosmological horizon in its interior, therefore including both a local and a 
cosmological regions (as opposed to considering only the local region) . This allows us to determine decays 
for initial data which lead to uniform exponential decay in time of the solution. 



1 This can be seen from the fact that operator T in equation ( 1131) , whose fixed points are the solutions of Q , is a constant 
operator for A = 0; when perturbing to the nonlinear problem this operator becomes a contraction for small initial data. 
2 In this case the operator T is not even a contraction in the full domain. 



Numerical evidence for pointwise exponential decay can be found in [5] , and references therein, where 
higher spherical harmonics are also studied, as well as the non-linear system. Yagdjian and Galstian [T2] 
constructed the fundamental solutions of (TTJ) in de Sitter spacetime and proved exponential decay of 
certain homogeneous Sobolev L p norms, 2 < p < oo. Also along these lines, Ringstrom [11] obtained 
exponential decay for non-linear perturbations of locally de Sitter cosmological models in the context of 
the Einstein-nonlinear scalar field system with a positive potential. Finally, exponential pointwise decay 
in the local region between the black hole and the cosmological horizons in a Schwarzschild-de-Sittcr 
spacetime follows from the papers by Dafermos and Rodnianski [51 E] (see also [T]). 

2 Christodoulou's framework for spherical waves 

Bondi coordinates [3j (u, r, 9, if) map the causal future of any point in de Sitter spacetime isometrically 
onto ([0, oo) x [0, oo) x S 2 ,g), where 



1- -r 2 ] dir ■ 2diulr + i J dif 



(2) 



with d£l 2 the round metric of the two-sphere (cf. Figure [TJ . 

i J+ 




Figure 1: Penrose diagram of de Sitter spacetime. The lines u = constant arc the outgoing null geodesies 
starting at r = 0. The point i corresponds to u = +oo, the cosmological horizon H to r — \ ■?■ and the 
future null infinity ,J? + to r = oo. 

In these coordinates the wave equation 



□ 3 = o ^ c\, (y-detfoja^j = o , 

for spherically symmetric functions, dg(f> = dtptfi — 0, reads 



dr \duj du \ 3 / dr 2 
Following Christodoulou [3] we consider the change of variable 

h ■= §~ r M ■ 



If we assume that 



it immediately follows that 



lim r<b = 



>ri=° 



(3) 



, i r d<j> dh h-h 

= h := - / h (u, s) as and — = — = 

or or r 



(4) 



Moreover, assuming that the crossed partial derivatives of r<p commute, we see that ([3]) is equivalent to 

Dh=--r(h-h), (5) 

where D is the differential operator given by 

du 2 V 3 J dr 

3 Main result: statement and proof 

Our main result is the following 

Theorem 1. Let A > 0. Given ho £ C k ([0, oo)), for some k > 1, the problem 



Dh= -±r(h-h) 
h(0,r) = ho(r) 



(6) 



has a unique solution h e C ([0,oo) x [0,oo)). 
Moreover, ?/||/io||c° is finitqj then 

\\h\\c<> = \\ho\\c° ■ (7) 

Also, if ||(1 + r) p d r ho\\co is finite for some < p < 4 and H < 2'i/4 then 

||(1 + rfe Hu d r h\\ c o < ||(1 + r) p 9 r ft ||cc , (8) 

and, consequently, there exists hEW. such that 

\h(u,r)-h\ < (l + r) n(p) e- Hu , (9) 

with 

f \ f , 2<p<4 
»(P) = | 2 f o< p <2 ■ ^ 

Remark 1. The powers of 1 + r obtained are far from optimal. Since we are mainly interested in 
understanding whether the decay in u obtained by this method is uniform in r, we were only careful in 
computing precise estimates for 2 < p < 4, which is enough to establish uniform decay for p > 2 (if p > 4 
the p = A result applies, and in fact it does not seem to be possible to obtain a stronger decay in r for 
d r h). For p < 2 our method does not provide uniform decay, but it is not clear if this is an artifact of 
these techniques or an intrinsic property of spherical linear waves in de Sitter. 

Proof. For h € C°([0, oo) x [0, oo)), we have rh £ C°([0, oo) x [0, oo)), and so we can define F(h) to be 
the solution to the linear equation 



D{T{h)) = -±r(T(h)-h) 
T(h)(Q,r) = ho(r) 

The integral lines of D (incoming light rays in de Sitter) , which satisfy 



(11) 



dr 1 /„ A ,\ ,„ . 

^=-2( 1 "3 r J' ^ 

are characteristics of the problem at hand. Integrating (fTTl) along snch characteristics we obtain 

T{h)(u 1 ,r 1 ) = h o (r(0))e-^f° lr(s)ds + - / r(v)h(v, r(«))e"» ^ r ^ ds dv , (13) 



where, to simplify the notation, we denote the solution to (1121) satisfying r(u\) — r% simply by s i— > r(s); 
we are dropping any explicit reference to the dependence on (u\, ri), but it should be noted, in particular, 
that r(0) is an analytic function of (iti, ri). 



^Recall that if / : X — ► R is continuous and bounded then ||/||co = sup^ex 1/(^)1- 



Given U,R>0, let C^ R denote the Banach space (c° ([0, U] x [0, R]) , || • || c o a J , where 



l/lk 



sup |/(«,r)| 



(•u,r)6[0,C/]x[0,-R] 



(14) 



Let r c := «/^ be the unique non- negative zero of 1 — ^r 2 (see ([2]) and <jT2j)). The non-decreasing 

behavior of the characteristics satisfying r\ > r c shows that the restriction of J 7 to C[} R is well defined 
for all R> r c . In fact: 



Lemma 1. Given U > and R>r c := y x, J 7 contracts in Cjj R . 



Proof. Fix [/ > and i? > r c . Then 



||J"(/ii)- J"(/i 2 )|| c p = sup |J"(/ii)(Mi,ri)- J"(/i 2 )(ui,r-i)| 



(ui,ri)6[0,C/]x[0,H] 

< sup (^ / 1 r(w)|/ii(w,r(w))-n 2 (w,r(w))|e-^^ lr(s)ds dw 



< 



sup 



A '"' ; 
111 A 



(«i,ri)e[0,C/]x[O,iJ] LJo 



< sup 
(«i,n)e[0,£/]x[0,fl] 

< sup 

(«i,ri)G[0,tf]x[0,ii] 



I/" 1 |r(«)e-* J'" 1 ^ ds d«} .\\hi- h 2 \\ c 

/[e-tC lr Wd«l Ul I 

I L Jti=oJ 



sup <— / I hi (u, s) — fi2(u, s) | ds 

(u,r)6[0,!7]x[0,il] I r JO 



{l-e-^^ r{8)d8 }-\\hi-h4c° UiR - 



Throughout, to obtain estimates, and in particular to estimate a, one needs to consider three (causally) 
separate regions, naturally corresponding to the bifurcations of (|12p : the local region (r < r c ), the 
cosmological horizon (r = r c ), and the cosmological region (r > r c ). However, since the computations 
are similar we will only present the details concerning the most delicate case, r > r c . 

The solution to (TT21) satisfying n = r{u{) > r c := a/j, is given by 



r(uWiooth(i^(c-«) 



(15) 



where 



c = u 1 +2W — arcoth V/t^i 



(in particular c > u\, and so (fTS")) is well defined for < u < ui). It follows that 



r(s)ds = 



fcoth (W|(c-a) U 



2— In 
as 



sinh (^(c-*) 



ds = In 



sinh 



1 /A 



f(c-«i)) 



sinh (Wf c 



and consequently 



sinh 2 U\l^{c -mj\ 



cosh 



(l^A (c _ Ml) ) 



sinh" 



V2y 37 



sinh 



(iy/Jc) coth 2 (|^/f(c-«i)) 



cosh(a(c — Ui)) 1 



sinh (ae) 2oti 

-2cmi 



a(o— «l) I g-Q(c-Ui) 



d qc /,-ac 



4a 2 r 2 



> 



4a 2 r 2 



where a := \\ ^r- Define 



(U,R):= sup (l - e~$ ^ r(s)ds ) 



(«i,n)e[0,£/]x(r c ,il] 



< sup i___ T . < 1- T — 52— <1. 

(«i,ri)e[0,C/]x(r c ,fl] V 4a r l / \ A R<> J 



Similar computations give 

sup (l - e^fo 1 -(.)-) < f 1 - U^L ) < 1 , 

(ui,ri)£[0,l/]x[0,r e ) V V ' 

for the local region, and 



a hor := sup (l-e-^o" 1 ^< 1 _ e -Vf^<i ) 

«i£[0,C/] V 7 

along the cosmological horizon. Finally er = max{o7 oc , cr^or, o- cosm \ < 1, and the statement of the lemma 
follows. 

□ 

By the contraction mapping theorem [5], given U > and R > r Cl there exits a unique fixed point 
hu,R € C[} fl of J 7 . Uniqueness guarantees that in the intersection of two rectangles [0, U\] x [0, i?i]n[0, U2] x 
[0, -R2] the corresponding hjj lt Ri and hu z ,R 2 coincide. Consequently, there exists a unique continuous map 
h : [0, 00) x [0, 00) — > R such that /i = J"(/i), i.e., 

/ l ( Wl ,r 1 ) = / lo ( r (0)) e -^^ lr(s)ds + - / r(v)h(v,r(v)) e -^ ^ n r{s ^ ds dv , (16) 

in [0,oo) x [0,oo). Continuity of h implies continuity of rh, so we are allowed to differentiate (JTHJ) in 
the direction of D, which proves that h is in fact a (C ) solution of ((5|). Existence and uniqueness in 

C°([0,oo) x [0,oo)) follow. 

To see that a solution of ([5]) is as regular as its initial condition assume that ho G C k+1 , k > 0, and 
start by noticing that if h £ C k then rh and d r (rh) are also in C k . In particular for fc € C° we can 
differentiate (fTo]) with respect to iti to obtain 

^ = ^_( 6o(r („ ))e -*/.-<*.) + | (rS)K , ri) 



A [ Ul - d ( 

+ - / r(w)/i(u,r(w)) — (e 
3 J 9ui V 



_Aj-i r(s)ds 



This last expression shows that ft-o € C fc+1 and h <E C k implies J^- € C fc . The same reasoning works for 
the derivative with respect to n (note that although r(v) is given by different expressions according to 
whether n < r c , ri = r c or n > r c , it is still the solution of a smooth ODE satisfying r(u±) = r±, and as 
such depends smoothly on the data (m, 7*1)). Consequently, if h e C k+1 and h € C fe , then /i is in fact in 
C k+1 and the regularity statement follows by induction. 

To establish ([7]) first note that: 

Lemma 2. If \\h \\c° < Vo and \\h\\ c o < y , for some y > 0, then \\T(h)\\ c ° < yo- 

Proof. From (fTU| we see that 

\T(h)(ui, ri )\ < \\ho\\co e-*-To" lr W ds +||ft|| C o £ [ Ul r(v)e-$K lr ^ ds dv 

•J ,/n 



< J/0 (e~* ^ r W * + A p r(v)e -i ir r(*)«s dv \ = y Q . 



The last step follows by a direct computation, as before, or by noticing that since h = 1 is a solution 
to ([H]), with ho = 1, one has .F(l) = 1. 

□ 

Now consider the sequence 

h (u,r) = h (r) 
h n +l = F(h n ) 

We have already established that, for any U > and R > r c , h n converges in Cj}^. to h, the solution 
of ([5]) • Lemma [5] then tells us that 

IIMIc?, „ - IIMc < \\ho\\c° , and so \\h\\ c o = lim \\h n \\ c o < \\h \\ C o . 

U,H U,H n— ►oo U > H 

Since this holds for arbitrarily large U and R, the bound follows. 

We will now show that the estimate (JSJ) holds. First of all if h € C 1 we see that Dh and d r Dh are 
both continuous, and consequently Dd r h exists and its equal to d r Dh + [D, d r ]h\j. Using this last fact 
and equations ^j while differentiating ((3J) with respect to r we obtain an evolution equation for d r h: 

Dd r h = -2-r d r h . (18) 

Integrating the last equation along the (ingoing) characteristics, as before, yields 

d r h( Ul , n) = d r h (r )e- 2 -£f° U1 r « ds . (19) 

It is then clear that initial data controls the supremum norm of d r h. In fact, let 

d = \\(l + r) p d r h \\co. 
In the cosmological region (r > r c ), one has, after recalling (TTSl) , 

|(1 + r 1 ) J, e Hui drh(ui,ri) | = (l + r 1 ) p e J? " 1 9 r /jo(?-o)e-^- / 'o" lr(s)ds 

, ,' l + risinh(a(c-ui)) y gMi ^ sinh (q(c - m)) \ 4 ~ p 



1 + ro sinh (ac) / \ sinh (ac) 

where a = ^y ^ as before. Now, since c — u\ < c, then e - 2Q ( c - u i) > e~ 2ac , and 

sinh (a(c - m)) _ e «(c-«0 _ e -a(c-«i) 



oQC p-QC 



sinh (ac) e ac — e 

1 _ -2a(o-ui) (21) 



Also 



1 p— 2ac 



1 + n sinh (a(c - m)) 1 + ^ coth (a(c — ui)) sinh (a(c — Ui)) 
1 + r sinh (ac) 1 + ^ coth (ac) sinh (ac) 

sinh (ai(c — ui)) + ^ cosh (a(c — Ui)) 

sinh (a(c)) + ^ cosh (ac) (22) 

l+_to cosh (a(c- -»i)) 
~~ 77- cosh (ac) 

2a v ' 



< (2a + l)2e" 



QUI 



4 Here we are using the following generalized version of the Schwarz Lemma: if X and Y are two nonvanishing C 1 vector 
fields in R 2 and / is a C 1 function such that X ■ (Y ■ /) exists and is continuous then Y ■ (X ■ f) also exists and is equal to 
X-(Y-f)-[X,Y]-f. 



Therefore, if < p < 4 and H < 4a = 2^X73, we plug dHJ and ^ into (J2UJ) to obtain 



sup |(l + ri) p e ffMl 5 r /i(ui,ri)| < d sup 

(ui,n)£[0,U]x[r c ,R] («i,n)e[0,£/]x[r c ,fl] 



2 p (2a + l) p e^- 4 ") Ml 



<2 p {2a + l) p d . 



Similar, although simpler, computations yield 



sup |(l + n) p e ffui a r /i(ui,ri)| <16 sup |(1 + r x ) p d r h Q (r x )\ < 16d 

(ui,ri)e[0,t/]x[0,r c ] rie[0,r c ] 



(23) 



(24) 



for the local region. This proves © . 

To finish the proof of Theorem [T] all is left is to establish the uniform decay statement ([9]). Start with 

1 r 

\h(u,r) — h(u,r)\ < — / \h(u,r) — h(u,s)\ ds 



< 



<I 



r pr 



\d p h(u, p)\ dpds 



Js 

r t>r 



-H u 



Js 



e~ Hu 

(i + pY 



dpds 



< 



l+r > 2<p<4 

r-e _ff " , < p < 2 



These estimates for 2 < p < 4 are obtained by direct computation; they seem to be the optimal results 
which follow from this method. The remaining cases, with the exception of p = 0, are far from optimal. 
In fact, since we are mainly interested in a qualitative analysis, namely if the decay obtained is or not 
uniform in r (see Remark [T]), the results for p < 2 were obtained simply using ,.} , p < 1. 
Using ((5]) we then see that 



\d u h\ = 



< 



Dh + \(l- 


r 2 ) 


d r h 




~r(h-h) 


i 

+ 2 


(- 


y^ d r h 



< (l + r) n ^ e - Hu 



with n{p) as in the statement of the theorem. 

Now since d u h is integrable with respect to u, by the fundamental theorem of calculus, we see that 
there exists 



But 



lim h(u, r) = h(r) 



\h.(r 2 ) -h(ri)\ = lim \h(u,r 2 ) - h(u,ri)\ 



< lim 



\d r h(u, r)\dr 



< lim Ira —rile' 



-flu 



0, 



and, consequently, there exists h € K such that 



Finally 



h(r) = h . 

/■OO 

\h{u,r)-h\ < I \d v h(v,r)\dv 

■■ I (i + r ) n ^e- Hv dv < (1 + r ) n(p) e- Hv 



D 

Remark 2. TTie same calculation shows that given R > t/ie solutions of © satisfy \h(u,r)—h\ < e~ Hu 
uniformly for r € [0, i?], even i/ || (1 + r) p <9 r /io||c° is not finite. 



4 Boundedness and exponential uniform decay for spherical lin- 
ear waves in de Sitter 

We now translate part of the results in Theorem (TTJ) back into results concerning linear waves in de Sitter. 

Theorem 2. Let (M,g) be de Sitter spacetime with cosmological constant A and (u,r,0,tp) Bondi coor- 
dinates as in Section^ Let cf> — (f>(u, r) £ C 2 ([0, oo) x [0, oo)) be a solution^] to 



Then 



Moreover, if for some < p < 4 



n^ = o 



<sup|0 r (r^(O,r))| . (25) 

r>0 



sup 

r>0 



2 



(1 + r)"^ (/•"(()./■)) 



< oo , (26) 



dr 2 
then there exists <\> £ R such that, for H < 2-i/4 , 

\<f>{u,r) - <f>\ <(l + r) n ^e~ Hu , (27) 



where 



, . , , 2 <p<4 , . 

n(p) H 2 , 0<p<2 ' (28) 



Proof. Since <fr is a spherically symmetric C 2 solution of (fT]) we saw in Section [2] that h = d r (r(j>) satis- 
fies |J5]), with cf) = h. Applying Theorem ([T]) the results easily follow. □ 

Remark 3. Note that the bound (|25[) for <f>, unlike the bound ([7]) for h, depends on the derivative of the 
initial data, so we do have "loss of derivatives" in this case. 

Remark 4. Once more that the powers of 1 + r obtained are far from optimal, see Remark^ 

Remark 5. Lt should be emphasized that the boundedness and decay results are logically independent. In 
fact (|25p follows from (|7J| , which in turn is a consequence of a fortunate trick (see proof of Lemma 0) 
relying on the non-positivity of the factor of the zeroth order term in ([5]) (here, non-negativity of A) and 
the fact that T (|13[) is a contraction in appropriate function spaces; in some sense one is required to prove 
existence and uniqueness of (O in the process. That is no longer the case for obtaining ©, from which 
uniform decay of <fi follows. 
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